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A.ls  TRACT 


Experimental  heat  transfer  studies  wore  conducted  on  full  scale 
ribbons  of  flat  circular  type  parachutes  having  a  geometric  porosity  of 
20.5  percent.  A  pressurized  subsonic  wind  tunnel  served  as  the  flow 
facility.  The  average  approach  flow  velocity  was  134  feet  per  second. 
The  measurements  were  made  in  the  Reynolds  number  range  of  one  million 
to  ten  million,  where  the  Reynolds  number  is  based  on  the  ribbon  width 
of  2.1  inches  and  on  the  velocity  and  temperature  of  the  flow  in  the 
slots  between  the  ribbons.  The  pressure  ratio  applied  to  the  ribbon 
was  varied  from  1.4  to  26.3*  Using  a  transient  energy  balance,  local 
and  average  heat  transfer  data  were  obtained  arid  compared  with  .available: 
analyses. 
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NOMKNCLA'nmJK 


convective  he.aL  transfer  area ,  ft* 


speed  of  sound,  ft/soc 


specific  heat,  iitu/lh|u  -  °H 


width  of  parachute  ribbon,  "slat"  ,  it 


O  O 

convective  heat  transfer  cociTicient,  Hlu/hi—J'u  ~  H 
thermal  conductivity,  litu/hr-f  t-°R 


Mach  number 


M!  Nussclt  number  based  on  stagnation  conditions 
ko 


Nussclt  number'  based  on  wall  conditions 
kw 

absolute  static  pressure,  Ibg/in^ 

rate  of  heat  transfer  per  unit  area,  litu/hr-ft^ 


distance  along  ribbon  from  stagnation  point,  ft. 


=  Reynolds  number  based  on  slat  width  arid  sonic  conditions 
R "  in  the  slot 


V 

|JV 


slat  half  width,  ft 


temperatures,  °F  or  °ll 


recovery  temperature,  c  l1'  or  °R 


thickness,  ft 


external  flow  velocity  outside  boundary  layer,  ft/ sec 
undisturbed  velocity  approaching  model,  i't/r.oc 
density,  lbm/ft^ 


time,  hr 


duperocrtpls 


denotes  sonic  conditions  calculated  from  the  measured  stagnation 


conditions,  assuming  an  isentropic  expansion  to  sonic  velocity 


devera!  symbols  defined  in  the  text  and  until  Imi  ir  ■  !' 
included  hi  l, his  list. 


ill 


NOMENCLATUIIK  (con I'd) 


denotes  adiabatic  waJl  conditions,  q  -  0 
average  value 

conditions  based  on  slat  width  D 

denotes  conditions  downstream  oi'  ribbon 

denotes  undisturbed  conditions  upstream  oi  model 

denotes  conditions  at  the  downstre;im  stagnation  point  of  the 
ribbon 

rel'ors  to  ribbon  material 

denotes  stagnation  conditions,  also  condition  in  settling 
chiimber  before  wind  tunned  contraction 

denotes  conditions  at  the  stagnation  point 

denotes  local  conditions  at  the  surface  of  the  test  ribbon 

denotes  conditions  at  edge  of  boundary  layer 


I.  INTRODUCTION 


A  growing  interest  in  the  use  of  parachutes  for  aerospace-vehicle, 
recovery  operations  has  directed  attention  to  the  aerodynamic  heat 
transfer  problems  involved  when  high  speed  aerodynamic  decelerators 
are  employed.  A  variety  of  retardation  devices  have  been  proposed  which 
satisfy  such  dynamic  design  requirements  as  filling  time,  opening  shock, 
drag  establishment  and  control  and  stability.  However,  as  the  flight 
Mach  number  is  increased,  a  common  problem,  aerodynamic  heating,  is 
encountered.  At  the  present  time,  there  exists  neither  a  fully  predictive 

theory  nor  sufficient  experimental  data  to  permit  reliable  calculations 

\/ 

/\ 

of  the  heat  rates  involved  in  the  several  deceleration  schemes  (l). 

In  general,  aerodynamic  decelerators  must  possess  a  large  drag-to- 
weight  ratio  and  must  be  capable  of  being  stored  in  a  small  volume. 
Therefore,  most  proposed  drag  producing  devices  are  relatively  thin  and 
flexible.  In  addition,  dynamic  .stability  requires  venting  or  porosity. 

As  a  result,  most  retardation  techniques  employ  porous  surfaces  involving 
high  temperature  plastic,  glass,  or  wire  mesh  cloth  in  conjunction  with 
other  venting  techniques.  The  heat  vransfer  to  these  elements  of  mesh  is 
a  common  problem  in  parachute  performance  calculations.  Alternately, 
experimental  heat  transfer  data  must  be  applied  to  a  variety  of  situations. 

The  flow  field  associated  with  a  parachute  results  from  the  inter¬ 
action  of  the  large-scale  phenomena  enveloping  the  complete  body  and  the 
small-scale  processes  associated  with  the  local  flow  about  a  single  element 
in  the  mesh  material  itself.  The  large  scale  phenomena,  (parachute  shape 
and  size,  velocity,  altitude)  which  vary  with  the  particular  design,  are 


Numbers  enclosed  by  brackets  refer  to  references. 


1 


generally  involved  with  the  aerodynamics  of  the  problem  and  contribute 
end  or1  boundary  conu it ion a  (pro scare  ratio,  Reynolds  number)  to  the 
small- sc; ilo  phenomena . 

A .  Jar ac  Scale  Phenomena 

The  range  of  interests  in  the  present  study  consists  of  altitudes 
from  sea  level  up  to  200,000  i'cet  and  flight  Mach  numbers  ranging  from 
2  to  ),  A  concave  hemisphere  canopy  is  a  typical  shape  considered  for 
supersonic  operations  although  specification  of  the  exact  geometry  is 
not  necessary .  In  previous  work  (2),  it  was  foil,  that  a  good  starting 
point  for  study  of  the  large-scale  phenomena  was  the  simple  concave 
hemisphere.  This  provides  a  flow  amenable  to  analysis  and  experiment. 

The  experimental  results  of  Reference  (3)  indicate  that  at  angular  loca¬ 
tions  up  to  73  degrees  from  the  stagnation  point,  the  static  pressure 
on  the  upstream  surface  of  the  concave  hemisphere  is  equal,  to  the  total 
pressure.  This  means  that  there  is  practically  no  flow  near  and  parallel 
to  the  surface  of  the  hemisphere  except  near  the  edge.  Of  course,  in 
supersonic  flow,  a  normal  shock  would  stand  ahead  of  the  hemisphere  and 
this  would  produce  a  completely  subsonic  flow  regime  inside  of  the  concave 
hemisphere.  In  the  absence  of  any  chute  porosity,  the  pressure  on  the 
upstream  side  of  the  chute  is  closely  approximated  by  the  total  pressure 
behind  the  normal  shock.  In  determining  the  pressure  on  the  rearward 
(downstream)  face  of  the  chute,  we  are  led  to  a  base  pressure  phenomena 
described  by  H.  H.  Korst  (/+)• 

The  discussion  in  the  previous  paragraph  illustrates  a  tcclmique 
for  computing  the  pressures  and  honco  the  pressure  ratio  across  a  chute 
with  no  porosity.  The  analysis  for  the  base  pressure  problem  has  been 
extended  to  l.lu!  case  of  flow  through  the  chute  with  small  momentum  only. 


Bineo  this  last  condition  is  rarely  mot  in  the  ease  of  practical  parachute 
configurations,  the  pressure  ratio  across  the  chute  openings  must  bn 
found  from  experiment . 

B.  Una 11  Beale  Bhenomena 

A  high  temperature,  high  density  flow  is  produced  inside  the  canopy. 
This  flow  passes  through  the  individual  openings  of  the  mesh.  Heat  is 
transferred  from  the  hot  gases  to  the  mesh  elements.  Therefore,  it  is 
the  distribution  of  heat  flux  to  the  surfaces  of  a  mesh  element  which  is 
of  primary  importance.  The  element  is  bounded  on  each  aide  by  an  opening 
(slot)  through  which  the  oncoming  flow  passes.  The  typical  ribbon  (slat) 
may  be  considered  as  a  flat  plate  aligned  normal  to  the  flow.  The 
approaching  flow  may  or  may  not  separate  upstream  ol'  the  ribbon.  This 
process  has  not  been  clearly  defined.  The  flow  passing  through  the  slots 
exhibits  properties  which  depend  upon  the  overall  applied  pressure  ratio. 

In  most  flight  applications,  the  applied  pressure  ratio  will  be  such 
that  the  flow  in  the  slot  may  bo  considered  as  a  sonic  flow.  Additional 
Cinergy  is  available  for  further  expansion  of  the  stream  as  it  emerges 
from  the  slot  ouch  that  localized  regions  of  supersonic  flow  will,  occur 
downstream  of  the  slot.  Between  the  supersonic  jets  emerging  froin  the 
slot  and  the  rearward  facing  wall  of  the  ribbon,  two  regions  are  found. 

The  first  is  a  conventional  free- shear  layer  in  which  the  slrcamwiso 
component  of  velocity  diminishes  from  the  value  found  in  the  jot  to  a  near- 
zero  condition.  Adjacent  to  the  shear  layer  there  exists  a  small  region 
of  reverse  flow.  There  is  a  circulating  vortex  in  this  region.  The  local 
pressure  in  the  separation  bubble  is  determined  by  tile  local  base  pressure 
phenomena  (the  flow  issuing  from  the  slot).  The  mean  near-wake  pressure 
which  exists  further  downstream  is  determined  by  the  gross  geometry  arid 

These  two  pressures 


the  frees  tre;un  Mach  number. 


are  probably  not  cijual. 


Therefore,  further  reeompreanion  racists  downstream  of  tho  jots  issuing 
J'rom  the  individual  slots.  The  local  base  pressure  behind  an  individual 
ribbon  is  therefore  determined  by  the  slot  pressure  ratio  and  the  reeoiti- 
pression  process  in  the  shear*  layers . 

Tho  distribution  of  heat  flux  over  the  surface  of  a  single  ribbon  is 
of  primary  .importance.  A  model  which  applies  to  the  general  heat  exchange 
process  between  the  shear  layers  and  the  rearward  facing  wall  of'  the  ribbon 
has  been  proposed  by  Korsl  (4)  and  Chapman  (5).  Related  experiments  arc 
discussed  in  References  (6)  and  (?). 

Considering  a  single  ribbon  in  the  absence  ol'  adjacent  ribbons,  the 
heat  flux  increases  from  the  "apparent"  stagnation  point  value  (at  the 
center  of  the  ribbon)  to  higher  values  toward  the  outer  edge  of  the  ribbon, 
according  to  laminar  analysis.  Experimental  evidence  (tf)  covering  the; 
forward  face  of  such  a  ribbon  confirms  this  distribution  of  heat  flux,  and 
in  addition,  the  variation  of  heat  flux  with  pressure  level  (for  a  fixed 
velocity  gradient).  When  several  ribbons  are  grouped  together  in  a  manner 
simulating  a  ribbon  parachute,  tho  heat  flux  distribution  on  the  front- 
facing  surface  of  a  typical  ribbon  is  influenced  by  tho  presence  of  the 
adjacent  ribbons  (2).  The  heat  flux  also  increases  at  locations  away  from 
the  stagnation  point  but  at  a  rate  somewhat  faster  than  a  single  ribbon. 

The  increased  heat  flux  near  the  edge  of  the  ribbon  appears  to  be  due  to 
a  diminution  in  the  thickness  of  the  viscous  layer  adjacent  to  the  ribbon 
slot.  This  diminution  is  produced  by  the  interaction  of  the  slot  flow 
with  the  ribbon  boundary  layer. 

Recent  experiments  (6,7)  suggest  that  adjacent  to  the  separation 
bubble,  i l  laminar,  boundary- Jaycr-likc  flow  exists  along  the  rear  surface 
and  Ui  is  serves  as  the:  major  barrier  to  the  exchange  of  energy  between 

k 


the  free stream  and  the  surface.  In  the  case  of  a  flat  plate  aligned 
normal  to  the  flow,  this  thin  film  exhibits  the  properties  of  a 
stagnation- point  boundary  layer  of  the  laminar  type;.  The  heat  transfer 
processes  along  the  forward  and  rearward  faces  of  the  ribbon  appear  to 
be  governed  by  a  similar  mechanism. 

Experimentally  (7),  the  magnitude  of  the  heat  fluxes  to  the  forward 
face  are  found  to  be  considerably  larger  than  those  of  the  rearvrard  face. 
This  is  primarily  duo  to  the  significant  difference  in  pressure  level 
which  occurs  between  the  front  and  rear  surfaces.  The  specific  value  of 
the  ratio  of  heat  fluxes  on  the  forward  and  rearward  surfaces  varies  with 
the  freestream  unit  Reynolds  number  since  base  pressure  phenomena  depend 
on  the  character  of  the  shear  layer  (laminar,  transitional  or  turbulent). 

The  greatest  source  of  uncertainty  in  engineering  calculations  of 
the  energy  exchange  process  to  the  mesh  elements  lies  in  the  determination 
of  the  convective  heat  flux  as  characterized  by  a  heat  transfer  parameter 
Nu.  In  the  present  instance  the  heat  transfer  parameter  is  a  function  of 

Reynolds  number,  Re  }  and  the  pressure  ratio  across  the  ribbon,  \\/\'.} .  Thu 

■)(* 

inter-relationships  between  Nu,  Re  ,  and  1'^/lVj  are  sufficiently  complex 
that  they  must  be  determined  experimentally.  It  is  convenient  bo  define 
a  ReynoJ.dc  number  based  on  the  velocity  and  flow  properties  in  the  sonic 
orifice  and  the  slat  width 

in*  .  P"  a*  !>  (I) 

d  - — 

u 

In  the  flight  case  the  static  temperature  T  and  pressure  are 

determined  by  the  flight  altitude.  A  total  pressure  p  and  temperature 
T0  are  determined  once  the  flight  Mach  number  is  prescribed.  Since  the  how 


shock  is  assumed  to  be  normal,  Urn  bo  ha  I  (ircuuui'o  behind  the  normal  shock 
P!  (inside  o.f  the  canopy)  in  determined  from  p_.  and  M  .  1  The  total 
density,  behind  the  bow  shock,  is  computed  from  Tq  and  'the  proper 

equation  of  state.  Finally,  the  sonic  density  and  speed  are  determined 
from  an  isentropic  expanison  J'rom  id,  f  to  p"  and  T  .  The  an por script  * 
refers  to  the  slot  condition  where  a  sonic;  speed  exists  for  supercritical 
pressure  ratios.  The  viscosity  n"  is  a  known  function  of  the  temperature 


The  heat  transfer  parameter  Mu  is  based  on  the  stagnation  conditions 
of  the  approaching  flow  and  relates  the  wall  heat  flux  to  the  difference 
between  the  wall  temperature  Tw  and  the  adiabatic  wall  temperature  T  . 


m  -  hi.) 
Mu  r-  — 

ko 


q. 

T  ,,  -  T, 


Notice  that  k  is  the  thermal  conductivity  of  the  air  at  the  stagnation 
temperature . 

1 1 .  MX  I'KKIML'INTAI.  KQUltWiNT  AND  I’RUUKJJUHKl 

A .  I'Vu'.il  :i.ty 

The  experlmenta]  facility  used  in  the  present  experiments  appears  in 
the  photograph  of  the  experimental  setup.  Figure  2,  and  is  also  sketched 
in  Figure  j.  This  facility  is  located  at  the  University  of  Minnesota's 
Jtosemount  Id: search  Center.  I’rior  to  modification,  the-  wind  tunnel  was  a 
conventional  blowdown  supersonic  windtunnol  (v).  The  6  x  12- inch  supersonic 
vrindtunno.l  was  modified  by  removing  the  supersonic  nozzle  blocks  and  fabri¬ 
cating  new  0  x  8- inch  subsonic  test  section.  A  quick-acting  DeZ.urik  valve 


The  present  experimental  procedure  is  to  simulate  the  subsonjc  flow  cond.i- 
ti.ons  behind  the  bow  shock.  The  total  pressure  behind  the  normal  shock  j.;; 
properly  Simula. ted  by  the  wi.ndtunneJ  stagnation  pressure  (see  I.ion  1IA), 


was  installed  in  the  6- inch  inlet  pipe.  The  opening  time  of  this  valve 
is  slightly  less  than  one  second ,  which  is  of  importance  when  transient 
heat  transfer  measurements  are  made.  A  series  of  five  perforated  plates 
were  placed  at  11- inch  intervals  inside  the  18-inch  diameter  stilling 
chamber  and  a  16-mesh  screen  was  inserted  between  the  flanges  at  the 
downstream  end  of  this  chamber  to  even  out  the  flow.  The  plates  and 
screen  had  a  porosity  of  about  54  percent.  .Tests  reported  by  Baines  and 
Peterson  (10)  indicated  that  this  is  an  optimum  porosity  for  flattening 
out  non-uniform  velocity  profiles. 

A  short  round- to-rectangular  contraction  with  an  area  ratio  of  3.2 
followed  the  stilling  chamber,  and  the  supersonic  nozzle  blocks  normally 
were  installed  immediately  downstream.  These  blocks  were  replaced  by  a 
pair  of  two-dimensional  subsonic  contraction  blocks  with  a  contraction 
ratio  of  1.5.  The  contraction  contours  (cubical  arcs)  were  determined 
from  the  data  given  by  Rouse  and  Hasson  (ll)  on  a  cavitation  prevention 
basis.  The  second  contraction  provided  an  additional  acceleration  of 
the  flow  to  134  feet  per  second.  It  also  served  the  purpose  of  reducing 
the  nozzle  boundary  layer  thickness. 

The  compressed  air  system  consists  of  two  high-pressure  compressors, 
air  storage  tanks,  a  dryer,  and  associated  valving  and  piping.  One  com¬ 
pressor  is  a  two- stage  unit,  rated  at  250  cfm  at  600  psig  and  powered  by 
a  100  h.p.  electric  motor.  The  other  is  a  three- stage  compressor,  rated 
at  195  cfm  at  1500  psig,  and  driven  by  a  50  h.p.  electric  motor.  The 
compressors  pressurize  a  storage  volume  of  1460  cubic  feet  to  250  psia 
(18  atmospheres).  The  air  is  dried  to  a  dew  point  of  -40°F  prior  to 
storage . 
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Downstream  of  the  test  section  the  air  may  be  exhausted  either  to 
the  atmosphere  or  to  multiple  vacuum  tanks  whose  total  capacity  is 
22,750  cubic  feet.  The  vacuum  system  was  not  used  for  these  studies. 

Only  atmospheric  exhaust  was  employed. 

B.  The  Experimental  Model 

The  test  model  was  mounted  two  nozzle  heights  downstream  of  the 
nozzle  tangency  point  (see  Figure  3).  The  model  consisted  of  three 
slats  -  each  2.1- inches  high,  6- inches  wide  and  0.3 12- inch  thick  (see 
Figures  4  and  5 ) •  The  slats  were  mounted  in  such  a  way  that  they 
resemble  the  array  of  ribbons  in  a  ribbon  parachute.  A  single  ribbon 
spacing  of  0.537- inch  was  used.  This  applies  to  the  two  center  slots. 

The  outer  two  slots  were  0.337-inch.  This  value  was  obtained  by  selecting 
a  half-slot  height  plus  a  nozzle  displacement  thickness  correction.  A 
single  geometric  porosity  (open-to-total  area)  of  20.5  percent  was  examined”. 

This  porosity  was  selected  since  it  closely-  matched  the  design  throat  size 
of  the  tunnel  when  operated  supersonically;  i.e.,  design  subsonic  and  super¬ 
sonic  mass  flows  were  identical. 

At  the  maximum  upstream  pressure,  the  aerodynamic  loading  was 
approximately  200  pounds  per  square  inch  -  or  a  total  loading  per  slat  of 
approximately  one  ton.  The  test  slats  had  to  be  quite  husky  to  withstand 
this  loading.  The  three  slats  were  cantilevered  from  aluminum  "windows" 
on  each  end  as  may  be  observed  in  Figure  i+.  Optical  windows  were  embedded 
in  the  aluminum  windows  for  shadowgraph/s chlieren  observations  of  the  slot 
and  shear  layer  flows.  Structural  and  packaging  problems  limited  the  field 


Geometric  porosity  was  systematically  varied  in  Reference  (2)  from  20 
percent  to  84  percent.  Porosity  was  demonstrated  to  be  an  unimportant 
parameter  in  determining  the  heat  transfer  processes  in  a  subsonic 
ribbon  parachute.  This  conclusion  is  undoubtedly  still  true  in  super¬ 
sonic  flow  as  long  as  the  shock  remains  ahead  of  the  entire  parachute 
entrance . 
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of  view  to  a  4-inch  diameter  circle  located  slightly  downstream  of  the 
center  slat.  Two  views  of  the  ribbon  parachute  model,  as  installed  in 
the  wind  tunnel,  are  shown  as  Figures  l+a.  and  4b.  Certain  instrumenta¬ 
tion  details  are  also  visible  in  these  Figures  although  the  reader  is 
referred  to  Figure  5  for  additional  construction  and  instrumentation 
details . 

For  structural  reasons,  only  one  side  of  the  ribbon  was  instrumented. 
The  ribbon  orientation  was  fully  reversible.  The  basic  structure  of  the 
ribbon  was  a  solid  bar  of  type  304  stainless  steel.  Chambers  or  compart¬ 
ments  were  milled  in  the  bar  to  accept  the  instrumentation  leads  and  to 
provide  an  air  gap  element  of  insulation  between  the  forward  and  rearward¬ 
facing  surfaces  of  the  ribbons.  A  removable  constant-thickness,  stainless- 
steel,  instrumentation-plate  was  screwed  and  potted  to  the  main  load 
bearing  structure. 

Eleven,  36-gauge,  calibrated  iron-constantan  thermocouples  were 
cemented  at  regular  intervals  in  small  holes  drilled  in  the  internal  face 
of  the  instrumentation  plate  (see  Figure  5).  A  minimum  quantity  of  copper 
oxide  cement  was  used  in  the  installation.  Small  rectangular  grooves, 
0.025-inch  deep  and  0.040-inch  wide,  were  milled  in  the  stainless  steel 
instrumentation  plate.  These  grooves  were  filled  with  Resiweld  epoxy 
cement  and  the  epoxy  was  shaved  flush.  Two  additional  grooves,  visible 
in  Figure  4,  were  provided  on  the  exposed  face  of  the  instrumentation 
plate.  These  grooves  blocked  the  conductive  heat  flux  by  reducing  the 
effective  thermal  conductivity  of  the  surface  material. 

Two  small  tab  elements,  which  covered  a  portion  of  the  slot  region, 
were  brazed  to  the  edge  of  the  main  instrumentation  plate.  The  tabs  and 
the  main  plate  were  of  identical  material  and  of  equal  thickness.  A  single 
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thermocouple  was  installed  in  the  center  of  the  rear  face  of  each  tab. 
Hollow  pressure  chambers  were  soldered  to  the  sheltered  face  of  the 
instrumentation  plate.  The  pressure  taps  were  drilled  (d  =  0.020- inch) 
through  the  plate  into  the  chambers.  A  typical  pressure  tap  installation 
is  shown  in  Figure  5 • 

III.  EXPERIMENTAL  STUDIES 
A.  Uniformity  of  Approach  Flow 

As  the  airflow  enters  the  wind  tunnel  stagnation  chamber,  It  under¬ 
goes  a  sharp  90  degree  turn  (see  Figures  2  and  3)*  The  perforated  plates 
and  the  screen  were  installed  to  remove  any  resulting  nonuniformities. 

A  15-tube  velocity  survey  rake,  Figure  7,  was  fabricated  to  sense  the 
uniformity  of  the  approaching  flow  velocity.  The  test  ribbon  assembly 

X 

X 

was  moved  downstream  to  the  last  window  location'.'  The  velocity  survey 
rake  was  installed  in  the  abandoned  test  plane.  Runs  were  carried  out 
at  upstream  pressures  of  104  and  74  psia.  These  pressures  exceeded  the 
limits  of  available  manometric  instrumentation.  A  special  manometer, 
visible  in  Figure  2,  was  constructed  using  high  pressure,  translucent 
plastic  tubing.  Water  colored  with  food  dye  was  used  as  manometer  fluid. 
The  specific  gravity  of  this  combination,  as  measured  by  a  hydrometer, 
was  0.999  at  77°F.  A  great  deal  of  difficulty  was  encountered  in 
accurately  measuring  the  differences  in  impact  pressure.  For  example, 
at  a  stagnation  pressure  of  74  psia,  a  typical  difference  in  impact 
pressures  was  l/4-inch  of  water  (0.01  psia).  The  combined  effects  of 
high  pressure  level  (74  psia),  low  pressure  difference  (0.01  psia),  short 


This  was  necessary  since  the  test  ribbons  provide  the  choking  mechanism 
required  to  establish  a  fixed  value  of  the  approach  velocity. 
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wind  tunnel  running  time  (30  to  50  seconds),  and  the  various  response 
rates  of  the  individual  tubes  of  the  rake,  were  such  that,  initially, 
the  manometer  board  was  emptied  of  fluid  on  nearly  every  run.  The 
response  problem  was  examined  in  bench  tests  in  whi.ch  the  rake-manometer 
combination  was  connected  to  a  pressure  manifold  to  which  a  100  psia  test 
pressure  could  be  suddenly  applied.  Satisfactory  pressure  readings  were 
obtained  by  measuring  the  pressure  difference s  ' between  tubes  of  nearly 
equal  response  rates. 

The  approach  flow  velocity  distribution  data  are  summarized  in 
Figure  8.  All  velocities  are  within  one  percent  of  the  centerline  velocity 
with  the  exception  of  the  four  points  taken  near  the  nozzle  walls.  These 
latter  points  all  yield  low  indicated  velocities  which  are  undoubtedly  due 

X 

to  the  nozzle  boundary  layers.  The  average  approach  velocity  was  found  to 
be  134  feet  per  second.  The  uniformity  of  the  approach  flow  seemed  to  be 
within  tolerable  limits. 

B.  Ribbon  Pressure  Ratios 

The  pressure  ratios  applied  across  the  ribbon  elements  are  presented 
in  Figures  9  and  10  as  a  function  of  the  slot  sonic  Reynolds  number  Re^. 
This  quantity  is  defined  as  (12 ) 


Re 


!D 


p~'ca~'  D 


•k  p°  u 

372" 


D 


or 


2.270T"  x  10 
T*  +  198.6 


-8 


Re 


D 


-  5-U9  x  10  (psia) 


(3) 


(4) 


for  a  stag-.ai.ion  temperature  of  70°F  and  a  ribbon  width  of  2.1- inches. 

The  operating  range  of  the  wind  tunnel,  as  determined  by  a  useful  running 
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20  -  140  psia.  Two  downstream  pressures  were 


tamo  oJ'  30  seconds,  is  P, 


measured  for  reference.  The  first,  is  the  pressure  measured  at  the 

downstream  stagnation  point  of  the  ribbon.  The  second,  was  measured 

using  a  wall  static  pressure  tap  installed  in  the  centerline  of  the  last 

window  i  oii.'i  t.  inn  (see  Figure  3 ) .  p  a  ^  should  correspond  to  the  pressure 

p  measured  properly  on  an  actual  parachute.  This  occurs  several  ribbon 

widths  downstream  of  a  ribbon  element. 

bince  the  wind  tunnel  discharges  into  the  atmosphere,  at  low  flow 

rates  P’ox^^  is  nearly  atmospheric.  A  sonic  or  critical  pressure  ratio 

is  approximately  2.  Only  the  first  two  points  on  the  left  side  of  Figure 

V  are  subcritical.  All  others  are  supercritical.  This  point  is  clear  in 

Figure  10.  A  comparison  of  Figures  6,  9,  and  10  is  quite  informative. 

The  critical  condition  occurs  between  photo  b  and  photo  o  of  Figure;  6. 

The  jets  from  the  two  slots  diverge  as  the  upstream  pressure  (and  the 

pressure  ratio)  is  increased.  The  two  diverging  jets  mo  t  at  an  Rc^  of 

approximately  3*5  million.  No  further  jet  expansion  is  possible  since 

the  jet  flow  completely  fills  the  channel  downstream  of  the  ribbons. 

Therefore  He  pattern  becomes  stable;  i.c.,  insensitive  to  further 

i  ncreases  in  the  upstream  pressure  I'.. .  The  pressure  ratio  J'/p 

x  1  exit 

remains  constant.  The  ribbon  porosity,  open-to- total  area,  is  20.5 
percent.  This  produces  a  geometric  area  ratio  of  4>B8,  iixpanding  ono- 
dimennionaiiy,  this  generates  a  downstream  Mach  number  of  3-15,  and  a 
3 ingnnt ion- Lo- static  pressure  ratio  of  46.  Actually,  the  shear  layers 
occupy  a.  portion  of  the  flow  area  downstream  of  the  ribbons.  The 
measured  pressure  ratio,  1  p/l g,  is  26,  which  yields  an  average  dowti- 
stre.-uii  da  eh  number  of  2. Yd. 


At  low  Ren  the  discrepancy  between  P  and  1’  . ,  in  a  base  flow 

UJ  2  exit. 

phenomenon.  Consider  the  vertical  dashed  lines  in  Figures  ')  and  10  at 
lie,  -  3-3  x  10°  (photo  e  of  Figure  6).  The  jot  flew  expands  to  a  Mach 
number  of  approximately  2.3  and  then  undergoes  a  shock  recomprocsion" 
to  I’exit.  According  to  Figure  9,  13,/P  --  2.8/12  -  0.23.  The  predicted 

value,  according  to  the  base  pressure  theory  of  Korot  (4),  is  0.24. 

At  large  values  of  Re^,  the  jets  completely  fill  the  duct,  the  flows 
are  parallel  and  exhibit  the  conventional  shock  diamond  pattern.  Apparently 
this  pattern  continues  for  some  distance  down  the  duct.  Beyond  Hop  —  U  x 
10^  t  j.  is  less  than  P„, 

*  X»  ,c 

One  important  conclusion  to  be  drawn  from  Figure  10  is  that  the  ?• 
ribbon  pressure  ratio  is  constant  at  a  value  of  approximately  26  (depending 
on  the  definition)  for  slot  sonic  Reynolds  numbers  greater  than  4  x  10^. 

This  result  applies  to  the  case  of  atmospheric  discharge.  In  another  sorioc 
of  tests,  P^^  was  increased  from  subatmosphoric  values  to  progressively 
larger  values  by  covering  the  exhaust  flange  with  stacked  layers  of  perfor¬ 
ated  plater,.  In  this  manner,  P^/P  it  was  varied  from  3  to  26  at  a  single 
Itejj  of  6.3  x  10^.  Those  latter  studies  will  be  discussed  in  Section  D.2  (b). 

Referring  again  to  Figures  9  arid  10,  it  is  necessary  to  point  out 
that  the  pressure  ratio  remains  fixed  for  Re”  >  4  x  1(/J ,  Any  heat  transfer 
trends  which  occur  at  Re^>4  x  1(/J  are  truly  due  to  Reynolds;  number  and 
not  due  to  pressure  ratio. 

The  phenomena  observed  in  Figures  6,  9,  and  10  will  occur  on  actual 

ribbon  parachutes.  The  Reynolds  numbers  at  which  the  several  events  woujd 

occur  arc  not  necessarily  related  to  the  present  values  of  Re" 

I)' 

The  shock  is  a  normal  shock  iri  photo  e  of  Figure  6,  This  quickly  revert.;; 
to  the  weak  shock  system  shown  in  photo  f . 


C .  Ribbon  Ourl'ace  Pressure  Distribution 

The  front  surface  and  rearward  surface  wall  static  pressure 
distributions  were  obtained  using  the  static  pressure  taps  ‘visible  in 
Figure  5 •  The  results,  plotted  in  Figure  11,  present  a  local  pressure 
at  location  x  divided  by  the  stagnation  pressure,  at  location  x  ~  0. 

R  Is  one- half  of  the  slat  height.  Notice  that  the  pressures  were 
measured  along  two  vertical,  lines  -  alternate  points  belong  to  different 
spanwiso  stations.  The  data  was  taken  at  small  values  of  x/R  in  view  of 
the  anticipated  correlation  of  results  with  stagnation  flows. 

The  point  x.  =•  0  represents  the  geometric  centerline  of  the  ribbon. 
For  supercritical  pressure  ratios,  the  pressures  along  both  the  front 
and  back  surfaces  reach  maximum  values  near  the  geometric  centerline". 
Pressures  along  the  front  surface  arc  considerably  move  uniform  than 
those  along  the  roar  surface 

The  actual  location  of  the  stagnation  point  was  obtained  from  largo 
scale  plots  of  the  pressure  distributions.  This  data  is  summarised  in 
Figure  12.  The  apparent  stagnation  point  nearly  always  fell  below  the 
geometric  centerline  for  both  the  upstream  and  downstream  surfaces.  The 
average  displacement  (U.06~ineh)  is  not  large  when  one  considers  the 
internal  diameter  of  a  pressure  tap  is  approximately  0. 03-inch.  It  was 
not  possible  to  move  the  apparent  stagnation  point  closer  to  the  geometric 

The  single  exception  occurs  along  the  rear  face  at  a  subcritical  pressure 
ratio.  The  application  of  a  stagnation  point  model  to  the  rear  surf. -i.ee 
doer;  not  apply  at  subcritical  pressure  ratios. 

The  ribbons  were  of  finite  thickness.  The  static,  pressure  taps  installed 
in  tin:  center  of  the  slot  width  indicated  local  Mach  numbers  wel I  in 
excess  of  one.  This  implies  that  the  sonic  line  lies  near  the  upstream 
face  of  the  ribbon.  A  typical  value  is  p  /p  0.22  when  IR  ■  i_G'j 

t/  ,  LUO  w  -L  *■ 

psi.a,  lie"  L‘  8.5  x  10  . 
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centerline  by  a  rumple  rotation  of  the  ribbon-  Apparently  this  slight 
flow  abnormality  is  a  property  of  the  facility. 


Dy  assuming  that  the  measured  pressures  arc  constant  across  the 
thickness  of  the  boundary  layer,  and  that  the  flow  outside  of  the  boundary 
layer  expands  isentropically  from  p  .  to  p  ,  the  local  velocities-,  jus; l, 
outside  of  the  boundary  layer  may  be  determined .  ^his  data  is  presented 
in  Figure  13-  Of  primary  importance  in  this  figure  is  the  fact  that  both 
the  upstream  and  downstream  velocities  gijow  linearly  with  distance  from 
the  stagnation  point,  This  phenomena  is  characteristic  of  stagnation  flows 
and  is  a  familiar  result  for  the  upstream  surface.  Figure  13  presents  the 
strongest  evidence  for  considering  the  flow  adjacent  to  the  downstream 
surface  as  a  stagnation  flow  -  recalling  that  this  occurs  only  when  super¬ 
critical  pressure  ratios  are  involved. 

Since  the  velocities  outside  of  the  boundary  layer  grow  linearly  with 

distance,  the  velocity  gradient  may  be  computed.  These  results  are  presented 
■K- 

in  Figure  l/4  .  Two  values  are  given  at  each  Re£,  for  both  the  forward  and 

rearward  surfaces.  The  asymmetry  between  the  upper  and  lower  values  is 
undoubtedly  clue  to  the  stagnation  point  shifts  summarized  in  Figure  12 . 

The  upstream  surface  velocity  gradient  ;i.s  independent  of  Reynolds 
number  for  supercritical  pressure  ratios.  For  the  fully  developed  jot 
flows  (Re^w  4  x  10^*)  the  downstream  surface  velocity  gradient  decreases 
with  Reynolds  number.  Over  the  narrow  range  3  •  1>  x  10^  <  Re^<  S.'j  x  If/*, 


constant 


..-1/2  -1/2 

Uu|)  =  constant  (It?) 


(5) 


Actually  the  square  root  of  the:  velocity  gradient  is  plotted  since  thin 
quantity  is  of  importance  i.n  stagnation  poinl.  heal,  transfer  calculation:'.. 


Tin;  measured  upstream  surface  stagnation  point  velocity  gradients 


may  be  compared  with  analyticai  predictions.  For  flow  over  blunt  bodies, 
the  first  term  of  the  Taylor  series  expansion  for  the  surface  velocity 
is  the  linear  relation 

tig  Gx,  ( 6 ) 

where  C  is  the  local  velocity  gradient-external  to  the  boundary  layer. 

For  reference,  the  stagnation  point,  velocity  gradient  is  expressed  by 


c  -  LH®.  (7) 

R 

where  U&)  is  the  subsonic  velocity  approaching  the  body  anti  R  is  the  body 
radius,  j)  is  a  constant  depending  on  body  geometry.  The  following  table 
illustrates  a  few  numerical  values  of  p. 


Geometry  p  Reference 

Circular  Cylinder  2  11 

Sphere  3/2  13 

Using  an  approach  velocity  of  132*  feet  per  second,  a.  body  radius  of  1 . Op 

12 

feet,  and  a  p  of  2  yields 


a  value  which  is  22*  percent  larger  than  the  measured  value. 

Methods  for  predicting  the  imixiinum  recirculation  aono  velocity 
(Ufa  for  the  rear  face)  are  being  developed .  For  the  present  experimental 
conditions,  reference  (?)  given; 


2 


P.i 


CT  tan  u 


00 


Ub  is  the  mean  recirculation  velocity 

Uj  in  the  mean  jet  velocity  downstream  of  the  slot  at  the  maximum  jet 
width  ==  2000  ft/aec 


pb 

is  the 

mean  recirculation  zone  density 

pj 

cf- 

■— > 

O 

P  j  ij 

mean  jet  density  correspond! ng  to  Ui,_u.  =  D 

=  ^b_ 

T  0  i 

•  J 

*  h 

T°.j 

,Tj 

cr 

i:*tho 

jot  mix! ng  layer  spreading  factor,  20 

«  1 

KMji 

a 

is  the 

jot  expansion  hall' -angle,  tan  a  =  0.2 

Substituting  these  values  into  Equation  3 .  0  yields 
=  200  foot  per  second 

T(2)(2fiO)  ’  „  HI  nec-1/2' 

y  !'l7o;? 
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I).  Heat  Transfer  Measurements 
1.  Forw  ard  — h'ac inn  Quid' ac a 

( a )  Recovery  Tamporature  Distribution 

The  30-'.>0  second  wind  tunnel  running  time  did  not  permit  a 
standard  "run -to -recovery"  which  normally  requires  30  to  Ji'j  minuton. 
Extrapolation  techniques  had  to  bo  adopted.  Two  tunnol  runs  wore  mado 
to  find  tho  one  thermocouple)  of  uhe  seven  installed  in  the  stagnation 
temperature  rake  which  responded  closest  to  tho  average  of  tho  seven. 
Pour  runs  were  than  made  to  datormino  tho  distribution  of  roc ovary 
temperatures.  Kor  these  runs  the  model  was  allowed  to  respond  thermally 


arid  to  approach  the  ml  i  abatif  condition.  Thu  indicated  model  temperatures 


T.  and  tunnel  total  temperatures  T  were  measured  as  a  tunc  Lion  of  time, 
n.  o 

The  extrapolations  Lo  aero  hunt  transfer  wore  made  by  plotting  T^/T0  versus 
;]/i.  and  examining  the  curve  for  the  limit  as  t;  increases.  The  extrapolations 
never  involved  a  temperature  change  of  more:  than  2°F  from  the  last  measured 
value.  The  upstn •;im  surface  recovery  temperature  ratios  are  given  in  Figure 
\  ‘j .  At  the  stagnation  point,  x/it  --  0,  the  recovery  temperature  should  equal 
the  free  sire;  uu  total  temperature  at  these  low  temperature  levels.  This  is 
not  always  the  ease  (see  Figure  l‘j,  x/ll  -  0),  a  phenomena  undoubtedly  due 
to  thermal  conduction  in  the  surface  of  the  model.  A  one  percent  change  in 
Ti/’J’o  corresponds  to  a  i)./(  degree  Falironhoit  cliange  in  the  indicated  recovery 
temperature.  The  relatively  low  recovery  temperatures  in  the  slot  reflect 
the  large  velocities  that  exist  there.  It  is  important  to  note  that  the 
equilibrium  temperatures  in  the  slot  are  lower  than  at  the  stagnation  point. 
The  steady  state  operational  temperature  of  a  parachute  is  determined  by  the 
conditions  at  the  stagnation  point  and  not  at  the  slot1'. 

(b)  Heat  Transfer  Runs 

The  general  heat  transfer  experimental  technique  is  described 
in  detail  in  Appendix  A.  The  tooliniquo  consisted  of  precooling  the  model 
by  packing  dry  ice  around  the  ribbons  ;incl  waiting  for  10  to  1J>  minutes 
until  the  surface  temperature  was  -108  -  1°F.  The  tunnel  was  started  rapidly 


In  I  lection  2,  the  discussion  of  Figure  28  points  out  that  the  heat  transfer 
coefficients  in  the  slot  are  6  or  7  times  those  at  the  stagnation  point. 
There  Core,  transient  operational  temperatures  are  determined  by  conditions 
in  the  slot. 


(in  loss  than  1  second  after  the  initiation  oi'  the  valve  actuator  '" ) .  The 
dry  ice  was  blasted  off  by  the  air  stro;un.  'Die  model  heated  up  und<;r  the 
direct  action  of  both  .forced  convection  and  the  heat  capacity  of  the:  instru¬ 
mented  akin  of  the  model.  Heat  transfer  coefficients  were  calculated  for 
several  times  using  liquation  A-1Z,  of  Appendix  A.  The  heat  transfer 
coefficients  were  extrapolated  to  Kero  time  (a  uniform  wall  temperature 
condition)  ignoring  data  taken  in  the  J’irst  seconds.  Model  thermal 

symmetry  was  verified  such  that  only  one-half  of  the  central  ribbon  was 
extensively  studied.  The  stagnation  point  shifts,  observed  by  means  of 
the  surface  pressure  distributions,  were  also  clear  in  temperature  distrib¬ 
ution  plots  made  at  several  time  conditions.  The  render  is  again  referred 
to  Appendix  A  whore  a  sample  run  is  ilJustrated  and  the  various  steps  in 
the  data  reduction  process  arc  illustrated. 

The  dimensional  upstream  stagnation  point  heat  transfer  coefficients 
arc  presented  in  figure  16.  The  data  was  taken  for  Reynolds  numbers  as 
high  as  ten  times  those  previously  obtained  (.10 ^ ) .  The  heat  transfer 
coefficients  increase  with  nearly  the  one- half  power  of  tho  Reynolds  number, 
a  situation  typical  oi'  laminar  stagnation  point  boundary  layers.  The 
scatter  in  the  data  is  approximately  -  Id  percent.  The  heat  transfer  and 


Muring  the  period  after  the  upstream  plug  valve  begins,  to  open,  there  is 
a  rapid  compression  of  the  air  in  the  stagnation  chamber  which  lasts  about 
1  second.  The  indicated  stagnation  temperature,  as  indicated  by  seven 
shielded  ^6-gauge  iron-constantan  thermocouples,  rises  sharply.  After  an 
indicated  temperature  rise  of  nearly  1A)°1'',  the  temperature  dropped  to  within 
a  few  degrees  oi'  room  temperature  and  then  decreased  at  a  rate  depending  on 
the  rate  that  mass  was  removed  from  the  storage  tanks.  The  initial  transients 
(valve  opening  time,  compression  in  the  stagnation  chamber,  downstream 
pressure  decrease  to  the  desired  subaliiiosphor.i.o  level)  were  clearly  visible 
in  the  data  trances  and  wort:  completed  in  1-iid  seconds  (with  the  exception 
of  the  aforementioned  stagnation  temperature  drop). 


IV 


and  pressure  data  at  the  upstream  stagnation  point  are  combined  in  a 

dimensionless  group  in  Figure  17.  The  parameter 
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is  equal  to  0.49  (13)  for  a  laminar  two-dimensional,  stagnation  point  for 
a  wall- to- stream  temperature  ratio  of  0.7*  Since  the  data  was  extrapolated 
to  time  x  =  0,  the  dry  ice  temperature  was  used  in  evaluating  the  gas 
properties  k^.,  5  w»  and  pw.  The  density  and  velocity  gradients  were  obtained 
from  the  pressure  measurements,  Figures  9  and  14,  respectively. 

The  data  agrees  with  stagnation  point  predictions  for  low  Reynolds 

/  -JC.  /  V 

numbers  (1  x  10  ^  ReD  —  r  x  10°).  At  higher  Reynolds  numbers  the  data 
falls  20-30  percent  above  the  two-dimensional  laminar  stagnation  point 
boundary  layer  analysis  of  Reshotko  and  Cohen  (13).  Tne  applied  pressure 
ratio  increases  with  Re£  in  Figure  17-  The  flagged  points  in  Figure  17 
represent  data  taken  at  Re^  =  6.3  x  10°  at  several  applied  pressure  ratios 
(also  see  Figure  27).  As  mentioned  in  Figure  3,  for  these  tests  with 
increased  downstream  pressure,  a  stacked  set  of  perforated  plates  covered 
the  exhaust  opening.  The  downstream  pressure,  pexit,  was  varied  by  changing 
the  overlay  of  the  perforated  plates.  The  flagged  points  represent  forward 
surface  heat  transfer  measurements  taken  over  a  range  of  Pq/Pg^t  of  5  to  25. 

The  local  distribution  of  heat  transfer  coefficients  on  the  upstream 
side  is  presented  in  Figure  18.  For  easy  reference,  the  local  values  were 
non-dimensionalized  using  the  stagnation  point  values  given  in  Figure  16. 

This  was  done  in  view  of  the  overall  aim  of  correlating  all  results  with  the 
stagnation  point  heat  transfer  which  may  be  predicted  with  precision. 
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The  data  of  Figure  18  has  been  corrected  for  stagnation  point  shift  such 

that  x  represents  the  actual  distance  from  the  stagnation  point  and  R, 

the  distance  from  the  true  stagnation  point  to  the  slot.  Some  typical 

faired  curves  demonstrate  that  the  local  distributions  are  not  universal 

but  are  strongly  dependent  on  Reynolds  number.  Local  heat  transfer 

increases  rapidly  away  from  the  stagnation  point.  Local  heat  transfer 

coefficients  were  observed  which  were  eight  times  the  stagnation  point 

value.  Further  attempts  to  correlate  local  heat  transfer  coefficients 

with  stagnation  point  values  were  abandoned  at  this  point  because  of  the 

uncertainties  involved  in  the  distributions.  The  curve  labeled  Re"  =0.99 

D 

x  10^  in  Figure  18  represents  a  subcritical  pressure  ratio.  This  distribu¬ 
tion  agrees  well  with  comparable  subcritical  measurements  reported  in 
reference  (2).  Higher  Reynolds  number  distributions  generally  yield  higher 
local  values  relative  to  the  stagnation  point  level y/ 

The  average  upstream  surface  heat  transfer  coefficients  are  presented 
in  Figure  19*  This  data  was  obtained  by  numerical  integration  of  the  local 
distributions  illustrated  in  Figure  18.  Notice  that  the  average  heat 
transfer  coefficients  increase  linearly  with  Reynolds  number.  This  is  to 
be  compared  with  the  normal  laminar  square  root  dependence. 

In  Figure  20,  the  average  Nussell  numbers,  defined  as 


(10) 


are  plotted  for  the  upstream  side.  The  characteristic  length  in  both 
Nussell  and  Reynolds  numbers  is  the  width  of  the  ribbon.  The  character¬ 
istic  velocity  in  the  Reynolds  number  is  the  velocity  in  the  slot.  This 
enables  a  comparison  of  present  heat  transfer  data  with  previous  ribbon 


parachute  data  (2 )  and  data  obtained  on  txibe  bank;;  where  ;i.L  .i;;  customary 
to  base  the  Reynolds  number  on  the  tube  d.LimcLor  and  on  the  velocity  in 
the  narrowest  cross  sec  tion  (14). 

The  suberiL.i.ea.L  data  of  Reference  (2)  in  presented  in  Figure  20  to 
illustrate  the  overall  Reynolds  number  dependence''’'.  Tim  increase  of  the 
slope  of  the  c.ui'vo  with  higher  Reynolds  numbers ,  first,  observed  here  in 
Figure  17 ,  has  been  observed  by  other  investigators  on  circular  cylinders 
and  tube  bundles  (14,  lb). 

2 .  Rearward  Facing  Surface 

(a)  Recovery  Temperature  Distribution 

The  recovery  temperature  data  is  summarized  in  Figure  21.  The 
sharp  rise  in  recovery  temperature  at  Ro^aj  2,5  x  106  may  be  correlated  with 
the  establishment  of  a  completely  supersonic  jet  flow  downstream  of  the  slots 
(see  Figure  6).'  The  primary  variation  of  the  recovery  temperature  ratio, 
T.IW/T0  is  with  pressure  ratio  and  not  Reynolds  number .  For  Ro£»  R  x  10^ 
the  constancy  of  f.lw/T0  re  fleets  the  constancy  of  pressure  ratio  (see 
Figure  10 ). 


The  solid  symbols  in  Figuro  20  represent  data  taken  in  two  different  test 
facilities  and  using  different  experimental  techniques.  The  agreement 
between  the  throe  experiments  is  reassuring . 


At  low  subsonic  Reynolds  numbers  a  vortex  system  is  formed  in  the  wakes 
behind  obstacles .  The  phenomena  is  characterized  by  periodic  vortex 
shedding.  A  suppression  of  the  vortex  sheet  has  been  accomplished  by  a 
splitter  plate  which  increases  the  l  anporaturc;  recovery  factor  r  r  (T 
-  Tg  )/(T0  -  Tg)  from  0.1  to  0.8-0 ,<j  (16).  At  transonic  and 

supex'sonie  spends  where  the  vortex  sheet  disappears,  tho  recovery  factors 
arc  again  0.8-0. 9*  Shadowgraph  observations  ol’  the  super  sonic,  case  reveal 
a  a Lab  I e  vortex  pattern  with  no  shedding.  Tho  expansion- shock  pattern 
stab!  I  isos  tho  flow  pattern  according  to  the  laws,  ol'  forbidden  signals. 


The  heat  transfer  coefficients  measured  at  the  downstream 


stagnation  point  <••  re  presented  in  Figure  22.  A  comparison  of  Figure:; 

16  and  22  reveals  that,  the;  downstream  stagnation  point  heat  transfer 
rates  exceed  the  upstream  stagnation  point  values  for  Rej.  <  6  x  1U(j. 

This  is  rather  surprising  when  one  considers  the  downstream  density  is 
considerably  smaller.  The  initial  decrease  in  h...  ,  with  Rc^  i.s  probab.ly 

due  to  the  density  decrease  (Figure  y).  The  sharp  rise  in  the  range  2  x 
1C6  4r  Re” 46x  iU6  might  possibly  result  from  the  increased  recirculation 
velocities  (7 )  >  The  fraction  ol'  the  shear  layer  flow  which  is  reversed 
decreases  slightly  as  the  pressure  ratio  (and  the  jot  Mach  number)  is 
increased.  .Since  the  density  in  the  recirculation  zone  decreases,  increased 
recirculation  zone  velocities  are  required,  liquation  _j . 'j  relates  a  stagna¬ 
tion  point  velocity  gradient  to  an  approach  velocity.  On  the  other  hand, 
the  velocity  gradient  measurements  on  the  rearward  facing  surface  (see 
Figure  14)  exhibit  a  decrease  in  magnitude  as  Re^  is  increased.  Therefore, 
both  the  density  and  the  velocity  gradient  are  decreasing  with  increasing 
Re'j'  while  the  heat  transfer  coefficient  is  increasing.  Those  behavior 
patterns  are  not  typical  of  ordinary  stagnation  flows. 

The  apparently  sharp  discontinuity  In  the  results  at  lie" 

6  x  10 1  is  not  understood  at  the  present  time.  No  discontinuous  process 
is  observable  in  any  of  the  photographs  or  graphs  at  this  point .  The 
downstre;un  stagnation  point  heat  transfer  data,  arc:  compared  with  two- 
dimensional  laminar  stagnation  point  analyses  in  Figure  2;3 .  The  velocity 
gradient  and  density  wore  obtained  from  Figure:'.  l)  and  1 4.  Additional 
experimental  resu'ts  are  presented  for  the  condition  of  elevated  exhaust 


pressure:'..  The  results  arc  not  in  oonip’l etc  agreement;  either  with  each 
other  or  with  laminar  stagnation  point,  analysis.  A  few  of  the  point:; 
compare  well  with  analysis  while  the  rent  do  not.  Two  possible  explana- 
Liiiiu;  are  trie  increased  turbulence  level  of  the  recirculating  flow  at  high 
Re"  and/or  that  the  recirculating  flow  pulsates  between  several  patterns, . 

The  distributions  of  heat  transfer  coefficients  on  the  downstream 
side  arc  given  in  Figure  2 U-  The  average  or  integrated  valuer,  arc  presented 
in  Figure  2i>.  The  combination  of  effects  produce  an  average;  heat  transfer 
coefficient  which  is  independent  of  Reynolds  number.  The  abnormalities  in 
the  stagnation  point  heat  transfer  variations  with  Roe  aro  masked.  This 
fact,  coupled  with  the  complex  distributions  given  in  Figure  24,  make  it 
difficult  to  compare  the  behavior  of  the  average  and  stag, nation  point  heat 
transfer  coefficients.  A  com,  nrison  of  Figures  ly  and  2'j  reveals  that  the 
downstream  average  heat  transfer  coefficients  fall  well  below  the  corres¬ 
ponding  upstream  average  values.  This  point  is  pursued  in  Figu.  -  2b. 

Faired  curves  through  the  upstream  and  downstrc;un  average  heat  transfer 
coefficients  are  given  over  a  wide  Reynolds  number  range.  At  low  Reynolds 
numbers  the  he,  .u  trnnr'V.r  on  the  front  side  is  greater  than  that  on  the 
back  side.  With  increasing  Reynolds  number  the  heat  flux  to  the  back  side 
increases  more  rapidly  and  .for  Rejj  -  50,000  it  is  as  high  as  on  the  front 
side  -  and  even  higher-  as  Rejj  is  increased  further.  This  is  a.  behavior 
pattern  which  agrees  quantitatively  with  tubes  rind  Lube  bundles  in  cross 
flow  (16). 

The  present  data  continue  the  trends  described  above  for  subscril icnl 
pressure  .ratios".  However,  for  critical  pressure  ratios,  the  downstream 
heat  transfer  becomes  progressively  less  than  the  upstream.  If  the 
stagnation  point  analogy  is  valid,  the  ratio  of  upstream  to  downstream 

The  pressure  ratios  for  the  two  lowest  Reynolds  number  points  of  the 
pi-'  S'-nf  experiment  are  subcritlcal. 

2  A 


average  heat  transfer  flux  should  vary  directly  witn  the  square  roof  of 
the  density  ratio  or  approximately  with  (jP;i /p^ ) •  As  ;in  example  choose 
He"  u  B  x  10°.  From  Figure  10,  (Pq/p^ )  ~  5 .  From  Figures  19  and  25, 

hyp/hj0W1  -  3Vf/.lP-0~3.  Thus  the  proper-  pressure  ratio  trend  Is  observed  . 

The  problems  are  not  completely  understood,  however.  A  summary  of 
data  taken  at  lie"  —  6.3 ,x  10^,  with  several  exit  pressures,  is  given  in 
Figure  27.  The  upstream  pressure  is  held  constant  in  this  figure  while 
the  downstream  pressure  is  varied.  The  pressure  ratios  are  all  super¬ 
critical.  The  upstream  heat  transfer  is  sensibly  independent  of  the 
downstream  exit  pressure.  However,  the  downstream  stagnation  point  heat 
transfer  coefficient  decreases  as  the  exit  pressure  is  increased.  From 
Figure  30  it  can  be  concluded  that  the  variations  in  did  not  affect 


The  final  figure,  Figure  28,  presents  measurements  of  the  slot  heat 
transfer  coefficient  taken  with  the  instrumented  surface  facing  forward 
(bottom  slot)  and  rearward  (after  rotation  on  the  top  slot). 

The  data  displays  the  slight  asymmetry  found  in  all  measurements 


and  points  out  that  the  extreme  heat  flux  rates  occur  in  the  slot. 
Undoubtedly  the  slot  heating  rates  depend  on  the  ribbon  thickness  or 
thickness  ratio.  Figures  15  and  28  summarize  the  slot  heat  transfer 
behavior  as  follows:  (1)  the  heat  transfer  coefficients  in  the  slot  are 
the  largest  observed  on  the  ribbon;  and,  (2)  the  equilibrium  temperatures 
in  the  slot  are  the  lowest  observed  on  the  ribbon,  blot  heat  transfer  is 


of  major  importance  in  transient  problems. 


IV.  CONCLUSIONS 


Kxperimcnts  have  been  performed  on  models  representing  the  individual 
ribbon:;  of  a  ribbon  parachute.  The  recovery  factor  mid  beat  transfer 
coefficient  results  obtained  are  valid  in  the  continuum  flow  regime  at 
I'randti  number  values  of  approximately  ().’/•  These  results  extend  the  region 
of  available  experimental  heat  transfer  data  on  ribbon  parachutes  to  include 
l£c;  of  1,000  to  10,000,000.  The  average  measured  Nusscit  numbers  are 
presented  in  figure  26.  It  was  found  that  the  maximum  local  heat  transfer 
coefficient  exists  at  the  edges  of  the  ribbon.  Also  the  maximum  local 
values  on  both  the  upstream  and  downstream  sides  also  occurred  at  the  edges. 
The  central  regions  of  the  upstream  and  downstream!  sides  exhibited  behavior 
similar  to  laminar  stagnation  point,  boundary  layers.  The  relative  minimum 
heat  flux,  rates  occurred  at  the  stagnation  point  for  each  surface-flow 
orientation.  Attempts  to  correlate  the  gross  ribbon  heat  transfer  character¬ 
istics  with  the  upstream  stagnation  point  heat  transfer  predictions  wero 
only  partially  successful. 

Large  differences  in  local  heat  transfer  behavior  were  attributable 
to  the;  onset  of  sonic  flow  in  the  slots,  between  the  ribbons.  The  applied 
pressure  ratio  was  demonstrated  to  affect  the  ratio  of  average  upstream  to 
average  downstrc.'im  heat  transfer  coefficients. 
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FIGURE  I.  SKETCH  Of  RIB30N  PARACHUTE  FLC.V  PATTERN'S 


FIGURE  3.  SKETCH  OF  THE  TEST  SECTION  OF  THE  PRESSURIZED  SUBSONIC 
WIND  TUNNEL 
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FIGURE  0.  UNIFORMITY  OF  APPROACH  FLOW  VELOCITIES 
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FIGURE  ! 2.  THE  SHIFT  IN  THE  STAGNATION  POINT  LOCATION  AS  INDICATED  BY 
SURFACE  PRESSURE  VEASURE.VENTS 


FRONT  UPPER 


FIGURE  14.  STAGNATION  POINT  VELOCITY  GRADIENTS  CALCULATED  FROM 
SURFACE  STATIC  PRESSURE  MEASUREMENTS 


RECOVERY  TEMPERATURE  RATIOS 


IGURE  15.  RECOVERY  TEMPERATURE  RATIOS  FOR  THE  UPSTREAM.  STDE 


UPSTREA.V  STAGNATION’  POINT  HEAT  TRANSFER 


ATMOSPHERIC  EXHAUST 


CORRELATION  CF  UPSTREAV  STAGNATION  POINT  HEAT  TRANSFER  AN 
PRESSURE  DATA  V/ITH  LAMINAR  STAGNATION  POINT  BOUNDARY  LAYE 
ANALYSIS 
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FIGURE  18.  DISTRIBUTION  OF  IIP  AT  TRANSFER  COEFFICIENTS  ON  THE  UP STREAM 


UPSTREAM  SURFACE  AVERAGE 


FIGURE  !  9,  AVERAGE  HEAT  TRANSFER  CCEFF I  Cl  E.'.'TS  CK  THE  UPSTREAV  SIDE 


DOWNSTREAM  RECOVERY  TEMPERATURE  RATIOS 


RECOVERY  TEMPERATURE  RATIOS  FOR  THE  DOWNSTREA,1 


FIGURE  22.  DO’VNSTREA.V  STAGNATION  POINT  HEAT  TRANSFER  COEFFICIENTS 


ATMOSPHERIC  EXHAUST 


CORRELATION'  OF  DOWNST.REAV.  STAGNATION  POINT  HEAT  TRANSFER 
AND  PRESSURE  DATA  WITH  LAV'INAR  STAGNATION  POINT  BOUNDARY 
LAYER  ANALYSIS 


DOWNSTREAM  SURFACE  HEAT  TRANSFER  DISTRIBUTIONS 


BUTTON  OF  HEAT  TRANSFER  COEFFICIENTS  ON  THE  DOWNSTREAM 


DOWNSTREAM  SURFACE  HEAT  TRANSFER  DISTRIBUTIONS 


DOWNSTREAM  AVERAGE  HEAT  TRANSFER  COEFFICIENTS 


COEFFICIE 


DOWNSTREAM  SURFACE  AVERAGE 


HEAT  transfer  data  at  the  center  of  the  slot  of  a  ribbon 
PARACHUTE  OF  EXAGGERATED  THICKNESS 


APPENDIX  A 


1.  Transient  Energy  Balance 

The  transient  technique  relies  on  the  heat  capacity  of  the  model 
skin  for  the  experimental  determination  of  convective  heat  flux.  The 
general  heat  balance  on  a  small  element  of  surface  must  include 


qbi  -  the  heat  flux  to/from  the  boundary  layer 
qsb  -  the  heat  stored 

Q.COnd  ~  the  heat  transferred  within  the  skin  by 
conduction 

1rad  “  the  heat  exchange  due  to  radiation  to  the 

wind  tunnel  (n)  and  to  the  model  interior  (m) 


term 

I 

II 

III 
IV 


9-bl  is  the  term  we  wish  to  determine  for  comparison  with  analytical 

X 

predictions.  To  generate  a  "good"  experiment,  qbl  should  be  the  dominant 


Figure  29  Energy  Balance  Notation 
The  general  heat  balance  is  written 


qbldAw  +  (qcondinAdin-  qcondoutAd^t)  -  qraddAw  -  9stdAw  (ll) 

The  temperature  is  assumed  to  be  constant  across  the  thickness  of  the 
surface  -  a  condition  referred  to  as  the  "thin-wall"  assumption,  i.e., 
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heat  colidueLi.u  into  :ui  o.li.writ  Of  surface  per  unit.  tune  is 

).  -mil  IV.-i)  .Mil  I  1  r  (  A  -I  v  ,d  A.,  -I  d^  T...  , 


-!•:  A,  ^  :m<l  t-hat.  out.  i.s  •  k  (A{,  -I  l_d  ^ )  (^W  A  fJx). 

1  c-ui  .bt  dx  \Lx 


Thu:;,  the  net  ml.p  of  lic.it.  conducted  into  the  element  is 


dA  ,  u  T 
d  w 


‘W«W(/‘<Ut  -  *  dx  -  km  dx 


d/\  -■  J.  iix,  A(j  -  tl  ;  (t  0 . O'/k- inch )  (14) 

The  heat  storage  tom  may  bo  written 

q-d'-^  °PmdW^  U5) 

where  dW  is  the  element  of  mass  lying  beneath  the  surface  area  dJ^. 

For  a  Mu i form  wall  thickness 

dW  -■  (>m  1  dx  t  (16) 

It  is  possible  to  obtain  reasonable,  but  not  precise,  estimates  of  the 
radiation  heat  fluxes.  Radiation  losses  to  the  test  section  walls  and 
nozzle  blocks  must  be  considered .  Radiant  energy  impinging  on  the  rear 
face  of  the  test  surface  which  originates  inside  of  the  model  (i.e.,  the 
support  plate)  is  also  of  importance.  The  radiation  heat  exchange  between 
the  test  surface  and  the  wind  tunnel  nozzle  -  nozzle  box  combination  is 


given  by 


-  Vn  *w  <T»4  -  Tn4) 


The  surface  A^;  is  completely  surrounded  by  the  nozzle  and  is  flat.  It 
can  be  stated  immediately  that  l'w_n  1.  The  model  surface  area  is 
small  compared  to  the  mean  test  surface  -  nozzle  separation  di:  ,unce. 
The  fraction  of  reflected  radiation  which  returns  to  the  model  is  small 
enough  to  be  neglected.  The  interchange  factor, £  is  then 


€n~0.6,  £w_n^0.1 


The  massiveness  of  the  nozzle  blocks  and  the  inch  thick  test  section 


walls  allowed  the  assumption  that  the  nozzle  temperature  T  was  equal 
to  the  room  temperature  during  the  brief  20  second  runs. 

The  emissivity  of  the  model  surface  was  estimated  to  be  0.2i|.  which 
is  the  value  for  iron  at  100°!"  freshly  rubbed  with  emery  paper.  As  the 
test  surface  is  aerodynamically  heated  by  the  external  flow,  it  receives 
varying  amounts  of  radiant  energy  from  the  support  structure.  In  this 
case,  the  radiation  geometry  consists  of  /two  parallel,  closely- spaced 
flat  plates .  As  before 

*^w-m  ~  ''w-nAv  “  -O )  (18) 

Again  the  geometric  factor  F  is  unity.  The  interchange  factor  for  two 
parallel  walls  is 


£ 


w-m 


,  £ 
€ 


w-m 


0.2k 


(IV) 


Combining  all  elements  of  the  heat  balance  results  in 


I 


%1  = 


111 


dx.2 


b4)  - 


(20) 


In  view  of  Equation  20,  the  following  is  a  discussion  of  the 
experimental  technique : 

The  model  temperature  was  preset  by  surrounding  the  model  with  solid 
C0o(T  =;  -109.8ol'')>  prior  to  starting  the  tunnel.  These  runs  were  made 
during  a.  typical  Minnesota  summer  day  such  that  the  tunnel  temperature 
averaged  80°F.  Therefore,  the  initial  temperature  potential  (190°F)  was 
the  same  for  all  tests.  When  the  model  temperatures  wore  stable  and 


uniform,  the  wind  tunnel  was  started .  The  dry  ice  was  removed  by  the 


tunnel  starting  blast.  Model  temperatures  increased  with  time  as  a 


result  of  the  convective  heating  process.  Each  term  in  Equation  20  in 

time  dependent.  Initially  only  the  heat  storage  term  and  the  model- 

nozzle  radiation  terms  contribute.  As  time  proceeds,  the  other  terms 

increase  in  importance.  The  equilibrium  situation  occurs  when  term  II 

reduces  to  zero,  and  the  model  is  at  the  adiabatic  wall  (T  )  condition. 

aw 

The  accuracy  of  the  methods  depends  to  a  large  extent  on  the 
frequency  with  which  reliable  thermocouple  readings  are  recorded  near 
the  beginning  of  a  run.  Modification  M24  of  the  Dymec  DY-2010A  Data 
Acquisition  by a tern  provided  25  channel  input  continuously  monitored  at  a 
precise  T,ute  of  5  channels  per  second.  Taking  data  early  in  a  run  implies 
that  no  surface  heat  conduction  corrections  arc  necessary;  a  fact  which 
greatly  simplifies  data  reduction.  The  aerodynamic  heat  flux  is  then 
obtained  for  the  isothermal  surface  temperature  case,  which  is  the  most 
basic  and  easily  specified  case  and  the  only  case  for  which  most  heat 
transfer  analyses  apply. 

Neglecting  terms  IT  and  IV",  Equation  20  is  written 

%1  “  h(T.lw  -  Tw)  -  Pfn  t  £1™  (21) 

m  dr 

with  h  and  T  hold  constant,  and  at  c  -  o,  T  -  T(0),  Equation  21  yields 
aw  w 

the  particular  integral 


T 

a  v/ 


T 


aw 


T(0) 


hf 


p  c  t 
m  p 
MU 


(22) 


The  m.’ix-hiium  radiant  heat  flux  w;i.s  found  to  be  10  Utu/ft^-hr.  Since  the 
initial  temperature  was  lb()°E,  the  radiation  heat  transfer  coefficient  of 
l.O/.lVO  0.05,  was  at  most,  ].//,  percent  of  the  convective  heal  transfer 
coefficient  and  could  bo  safely  neglected. 


The  left  side  of  Equation  22  represents  the  ratio  of  the  temperature 
potential  remaining  after  x  seconds  to  the  initial  temperature  potential. 
For  small  values  of  time  (x  P^O),  Equation  23  be  written 

Tw  -  T(0)  =  hv  (23) 

Ta„  -  T(0) 

which  applies  quite  closely  to  the  case  of  constant  applied  heat  flux 
(%1  =  constant)  and  generates  a  constant  rate  of  change  of  wall 
temperature  with  time. 

Writing  Equation  22  for  two  times,  x^  and  x g  yields 


h  = 


pc  t 
m 


(t2-  *1) 


.Ine 


(24) 


Equation  24  was  used  in  the  actual  data  reduction  procedure  since  the 
quantities  Tw^  ,  Tw^  ,  and  the  time  difference  (x^  -  x^  ),  are  known 

somewhat  more  accurately  than  Tw,  T(0)  and  the  absolute  time  -  as 
required  by  Equation  22  .  The  indicated  heat  transfer  coefficients 
obtained  using  Equation  24  were  extrapolated  to  the  condition  of  constant 
surface  temperature,  (t  =  0),  by  plotting  h  versus  the  consumed  temperature 
potential  ratio  (Tave-  T(0)  )/  (T&w  -  T(0)  ).  This  ratio  is  equal  to  zero 
when  t  =  0. 


The  analysis  of  a  sample  run  is  presented  in  Figures  30  through  33' 
The  temperature  traces,  the  basic  raw  data,  are  given  in  Figure  30  .  The 
initial  time,  x  =  0,  is  found  by  inspection.  The  variation  of  heat  flux, 
qbi,  with  surface  location  is  directly  proportional  to  the  trace  slopes 

W 


since  the  product,  Pmcn  t,  is  the  same  for  all  positions 


"'The  specific  heat  does  vary  strongly  with  temperature.  In  Figure  31  heat 
capacity  data  for  pure  iron  is  presented  which  should  apply  closely  to  the 
present  material,  type  304  stainless  steel. 
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The  stagnation  point  heat  transfer  rates  are  seen  to  be  the  smallest  of 
the  group  while  the  slot  rates  are  the  largest.  Heat  transfer  coefficients 
obtained  by  repeated  application  of  Equation  24  are  given  in  Figure  32 
together  with  a  visual  representation  of  the  effects  of  extrapolating  to 
a  uniform  wall  temperature  condition.  Experimentally,  there  seems  to  be 
a  time-temperature  region  in  which  conduction  effects  are  small  enough  to 
be  ignored.  This  region  does  not  include  the  slot  where  large  gradients 
in  temperature  are  rapidly  established. 

Finally,  Figure  33  illustrates  the  integration  of  the  local  heat 
transfer  values  -  a  process  which  yields  the  one  average  heat  transfer 
coefficient  for  the  run. 
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extrapolation  of  heat  transfer  coefficients  to  a  uniform 
WALL  TEMPERATURE  CONDITION 
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